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Summary 


We have investigated procedures for dealing with a variety 
of different sources of noise and systematic errors that are 
commonly found in Hadamard transform and also in other multiplex 
systems. We were able to set up criteria for comparing instruments 
in the Fourier and Hadamard transform families through the use 
of a commonly applicable mathematical formalism. The results 
obtained in this study have led to the completion of one Ph.D. 
thesis by a Cornell University graduate student, Ming-Hing Tai, 
and we have also published four journal articles under this same 
contract. The most recent one of these which has been accepj:ed 
for publication in 1978 appears here as Appendix. It contains 
most of the results we have been able to derive, in collaboration 
with Dr. N, J. A. Sloane, on sources of noise and errors, and the 
correction procedures that they call for. 



1, Introduction 




Our initial proposal submitted to NASA Langley in 1975 
undertook to investigate a number of sources of noise and error 
correcting procedures characteristic of Hadamard transform optical 
systems. Specific questions to which we proposed to address 
ourselves involved the following four topics — as listed verbatim 
in the original proposal. 

The questions of main concern are these: 

1. A spectrum or picture can always be 'framed' within 

\ 

a region of constant known brightness. In spectroscopy, for 
example, the ends of a spectrum can be blocked off, so that m 
spectral elements on each end are known to have zero intensity. 

If N measurements are then made to determine the N-2m unknown 
spectral intensities, there is a redundancy that enters because 
of the 2m constraints on the spectrum's end elements. Under 
these conditions, what is the optimum choice of the number m that 
will permit recovery of the best spectrum? Is there a better 
choice of an end pattern than a zero level constant signal? 

2, Continuing in the same vein, we can ask to what extent 
the end elements can be used to help reduction of spectral 
noise due to noise spikes in the data. We mentioned above that 
a noise spike can be eliminated with good results and replaced 
by a mean data value. However, if the end elements are known 

to be zero intensity, it should be possible to eliminate spikes 
more judiciously by forcing the replaced readings to replace 
zeroes in the end positions of the spectrum. 
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3. Turning to errors in the mask construction, we need 
to know 

i) What is the effect of random errors in the width of 
open and closed elements? 

ii) What is the effect when the open eiements, say, are 
systematically larger than the closed ones? This is an effect 
that can be brought about by excessive 'undercutting' in the 
photo-etch process frequently used in mask production. 

iii) What errors are brought about when the steps in the 
mask motion systematically differ by a small amount from the size 
of the encoding element? In a typical application, the mask may 
have a dimension measuring some N elements. At the end of the 
pass, however, the mask has only been moved some (N - c) elements, 
where, say, e < 1/2, What is the cumulative effect of this kind 
of error? How can it best be compensated? In practice stepper 
motors are driven by an accurate screw, and the errors discussed 
here are errors inherent in the construction of the screw. 

4. In addition we also undertook "a study of a number of 
important limitations of this family of instruments that have not 
been attacked to date. Some of these correspond to related 
limitations found in Fourier transform instruments." 

2. Results 

We have been .able to obtain an improved understanding of all 
of these effects, and have in fact been able to obtain analytical 
insights into a larger variety of sources of noise and error than 
had been originally anticipated. 


We start with item 4 listed above. Hero two publications 
have appeared in print (see section 3 below for publications 
1 and 2) . These two papers show the relative performance of 
Fourier and Iladamard transform spectrometers operated under 
identical dotector-noise-limited conditions. T\';o main conclusions 
emerge from these papers. First, the singly-encoded Hadamard 
transform systems encode with an efficiency that is roughly /2 
times higher than the corresponding efficiency obtained with a 
Michelson interferometric spectrometer. Second, the Michelson 
interferometer is less efficient than presented in standard texts. 
This second statement is not of the greatest importance because 
there had been some papers in the literature that had derived 
the correct performance figures; unfortunately, however, these 
papers had been widely ignored. Our contribution on this score 
was primarily to focus attention on the omission. 

Among the topics listed under section 1, item 3., systematic 
means for dealing with mask defects of the type 3 i) and 3 iii) are 
treated e.xhaustively in publication 4 {see Appendix A and the 
list in section 3 below) . In that paper written with Dr. Neil 
J. A. Sloane of the Bell Laboratories, Murray Hill, N.J., we have 
discussed: effects of faulty mask alignment; differences between 

slit width and step size; excessive gaps between encoding and 
and blocking masks, distortions produced by continuously moving, 
rather than stepped, masks, and so on. 

Item 3, ii) is a complicated problem that turned out to 
have an analytical solution, described in a previously published 
paper, written before the grant award had actually been made. 
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"Errors in Hadamard Transform Spectroscopy or Imaging Caused by 
Imperfect Masks" M.H. Tai, M. Harwit, and N.J.A. Sloane, Appl.Opt. 

2678 (1975) . 

Errors listed under item 2 of section 1 also are treated in 
paper 4 (again, see Appendix) . We were able to treat the effect 
of noise spikes, random noise and sinusoidally shaped noise froms. 
From these other types of noise can be analyzed by means of Fourier 
techniques. The question of obtaining more data points than are 
strictly necessary is one whose mathematical treatment bears the 
name of 'theory of singular weighing designs'. Again, the detailed 
analysis is difficult to summarize, but the raathemotical treatment 
is given in Appendix A, and we provide results there on optimum 
mask patterns that provide the best signal to noise ratios. 

Finally, turning to item 1, the theory of singular designs 
also applies here, and again optimum mask patterns can be selected 
for given conditions. Improvements in signal to noise ratio can 
in fact be obtained by constraining some of the values of the 
spectrum's end elements to take on a constant known brightness. 

Here the precise value of the brightness enters not di'rectly into 
the signal-to-noise ratio. It affects system performance more 
immediately in the requirements on dynamic range. By choosing 
the frame brightness to equal the mean brightness of the entire 
spectrum or scene, the dynamic range required by the data handling 
system can be made small, and this facilitates data handling. 



3 . Publications 


1. Fourier and lladamard transform spectrometers: a 

limited comparison, Ming Iling Tai and Martin Harwit, 

Applied Optics 2664 (1976). 

2. Fourier and lladamard transform spectrometers; a limited 
comparison II, Martin Marwit and Ming Hing Tai, Applied 
Optics, U, 3071 (1977) . 

3. lladamard transform techniques in spectroscopy and imaging 
spectroscopy, M. Harwit, in "Transform Techniques in Chemistry" 
edited by Peter Griffiths, Plenum Press 1978. 

4. Errors in Hadamard transform optics, N. J. A. Sloane, 

Martin Harwit and Ming Ming Tai, to be published. Applied 
Optics, 1978 

5. "A practical Hadamard transform spectrometer for astro- 
nomical application," Ming Hing Tai, Cornell University Ph.D. 
thesis and Center for Radiophysics and Space Research Report, 
CRSR 655 (208 pages) . 

4. Ph.D. Thesis Work 

A graduate student, Mr. Ming Hing Tai, was partially 
supported by the present contract while working on his Ph.D. 
thesis. His thesis work dealt with the performance of Hadamard 
transform spectrometers and the astronomical results that could 
be obtained with them, particularly in planetary observations. 
The thesis appeared as a Center for Radiophysics and Space 
Research (Cornell University) report CRSR 655, listed above. 
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Part of this 208 page compendium involves procedures for dealing 
with noise and systematic errors. 

5. Personnel 

Prof. M. Uarwit, Principal Investigator 
Mr. M. H. Tai, Graduate Research Assistant 
Mr. G. Melnick/ Graduate Reserach Assistant 
Mr, F. Dain, undergraduate Reserach Assistant 
Mr. G. Stasavage^ Research Technician 

6. Budget Status 

No uncoimnitted funds, as of November 20, 1977. 
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ABSTRACT 

This paper analyzes the distortion in Hadamard transform optical instru- 
ments caused by moving masks, incorrect mask alignment, missing measure- 
ments. diffraction, etc,, and describes techniques for reducing or eliminating 
this distortion. In a great many cases the behavior of the instrument can be 
characterized by a single matrix equation of the form t? ” TWa -r e, where the 
components of t) are the measurements, T is a matrix characterizing the instru- 
ment, \V specifies the mask configurations, a is a vector containing the unk- 
nown spectral intensities, and the components of e are small random errors. 


* Prewnt adclress; Langley Research Cenier, NASA, Hampton, Va., 2j665. 
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1. Introduction 

Hadamard transform optical systems have been widely studied during the past few years. 
By encoding the light with properly designed masks these instruments benefit from the so- 
called multiplex advantatie and achieve an increased signai-to-noise ratio. A variety of multiplex- 
ing instruments have been described in the literature, designed to improve images, spectra or 
both.'-'° 

Any optical technique that leads to improved performance is likely to have its limitations, 
Some of these limitations may be inherent in the design of the system, while others may be 
caused by practical difficulties' with the apparatus. In either case once such a limitation is recog- 
nized ways of working with Or around it can be e,xamined. The purpose of this paper is to 
describe a number of sources of error encountered in Hadamard transform optical instruments, 
to give techniques for reducing or eliminating these errors where such techniques have already 
been worked out, and to suggest some possible ways of dealing with errors that have not yet 
been fully studied. 

We shall mainly describe errors occurring in singly encoded spectrometers, although most 
of the discussion will apply equally well to imagers. Similar errors occur in doubly encoded sys- 
tems but have not yet been extensively investigated, 

2. Description of Instrument and Spectrum: Simplest Case 

In this section we consider a multiplexing spectrometer in its simplest mode of operation. 
We give a mathematical description of the instrument under the assumption that no errors 
occur except random fluctuations in the detector readings. Subsequent sections will consider 
various departures from this ideal behavior. 

2.1 Description of Instrument 

Let us consider then a spectrometer with a single narrow entrance slit and with a mask in 
the exit focal plane. It is assumed that the entrance slit is sufficiently narrow so as to be 
always filled by the incoming radiation. The input to the spectrometer can therefore be com- 
pletely specified by giving the frequency distribution of the radiation. It is most convenient to 
express this as a function F{v), say, of the wavenumber v. We further assume that over the 
operating range in which we are interested (he optical system between the entrance slit and the 
exit mask behaves like a linear system,’'. In particular if the input consists of monochromatic 



light at wavenumber and with unit intensity, so that F(i>) « then the distribution 

of intensity along the exit focal plane can be written as 

—Vo) . 

(We assume that distances along the exit focal plane have been calibrated in terms of 
wavenumber.) Thus //(v-vo) is the instrument's spectral display for monochromatic input at 
wavenumber vo. //Iv— vo) is called the unptilic response or point spread function^' of the instru- 
ment. Because of aberrations and diffraction //(v— v^) can never be concentrated at a single 
point, no matter how narrow the entrance slit is. By superposition the distribution of intensity 
along the exit focal plane produced by an arbitrary input Flo) is given by 

ck» 

C7(v) J* f {vq)//(v — vo)«/vo , (1) 

In other words C7(v) is the spectral display along the exit focal plane if the input spectttim is 
FM. 

An ideal instrument, free of aberrations and dilTraction, would have impulse response 
- 6(v— vq) and (from U)) 

C7(v) « FM . 

Thus the linear spectral display Glu) in the exit focal plane would be identical to the spectral 
distribution of the input. But in any real instrument Gin) Fin). 

or course strictly speaking we should write C7 as a function C7(.v) of the position, .v, along 
the exit focal plane, rather than as a function of wavenumber. But the latter notation is simpler 
and more suggestive. 

Suppose a single exit slit is used, extending say from wavenumber vj to wavenumber v; 
along the exit focal plane. .The detector output is given by the integral of the exit plane distri- 
bution over the slit area; 


T?(v|,vj)** J* Gi,n)di> + e 


« ^ J 

- J Fil>i^)dV(^ J //(v— vg)(/v 


e 


oo 

“ f f(Vo).S(Vo:Vl, Vj)rfVg + e , 


( 2 ) 


where e is the error due to detector noise, and 5(vo;v|, v;) is defined by 

5{vo;vj, vj) " J" H{n—i>i))du (3) 

■'1 

Vi) is called the slit j'unctioiP of the instrument; it is the total intensity of radiation 
passing through a slit extending from vj to V 2 in the exit focal plane, when the input is mono- 
chromatic radiation at wavenumber vg. 

Figures 1(a) and Kb) show two examples of F/(v-vo). Figure Ka) corresponds to an 
input slit which is wide enough for diffraction to be ignored, while Fig. 1 (b) shows the opposite 
extreme when diffraction effects dominate and //(v — vg) - rmcKv — vr,). where 
sine .X - (sin .v)/.v; The main lobe of /f(v~vo) has width 2d\ in Fig. I (a), or 2d: in Fig. Kb). 
The corresponding slit functions 5(vo;i'i.v;) are shown in Figs. Kc), Kd), assuming the exit 
slit has width v; — Vi equal to 2./) or Idi respectively. For e.xample in Fig. Kc) when 
ifQ «■ (vi+v 2)/2 the exit slit extends from vg - di to vg + d, and gathers ail the radiation, while 


if J'o “ (i^t+i' 2)/2 ±d{ only half of the radiation passes through the slit. 

Equation (2) is the fundamental equation that relates the output of the detector to the 
input spectrum FM when the exit slit extends from to t' - 

2.2 Description of Spectrum 

It is customary when measuring a spectrum with a dispersion instrument to divide the 
spectrum into n intervals, estimate the intensity in each interval, and then join the estimates by 
straight line segments. The result is a piecewise linear curve as shown in Fig, 2, Any spectrum 
can be approximated as closely as we please by such a curve, if n is taken large enough. We 
shall therefore start ojf with the .assumption that the input spectrum F(v) is a piecewise linear 

curve as shown in Fig. 2: the spectrum is divided into n segments 

corresponding to positions of the slits in the mask, and a, is the inten- 
sity at the center of (We assume that the scales have been chosen so that 9', and 9, are 

related by an equation of the form 9', *• 9, + A9^ -h • • • , where A is a constant.) 

The goal is to determine the n unknowns ao.fli, . . . as accurately as possible, in 
order to have an unambiguous expression for F(t>) in the first and last segments we assume 
that Fiv) is periodic oi.uside the interval 0o^,„ as shown by the broken lines in Fig. 2. Thus we 
define a~\ - o„_i and tf„ -> oq as the intensities in the segments <}-[9) to the left of Mi, and 
9„9„+\ to the righ^ of 9„^^9„^ respectively. Analytically we are assuming that FM is given by 

F{v) ■" a, + (n,+| — n,)(v — {l,4.| + i.V2f? (4) 

for 9,+\ ~~ b ^ V ^ 9,+\ + b 
where 2b « (?^+i — 9j is the length of each segment. 

2.3 Operation as a Multiplc.xcd Spectrometer 

Let us consider what happens when the instrument is operated as a multiplexed spectrom- 
eter, by allowing light from several slits to fall simultaneously onto the detector during each 
measurement. Suppose there are /;slii positions in all: 9r)9^,9\9'^, . . . ,9„-\0„. The pattern of 
open and closed slits is specified by means of the configuration matrix iV =■ (iv,,), where 
w^j - I if the y-th slit is open during the /-th measurement and w,, =« 0 othersvise, for 
0 < ij < /f— 1. In statistical terminology IK specifies a, wei^hins design}^“^^ The mechanical 
design of the instrument is simplified (as will appear in the next section) if IK is either a (right) 
circulant matrix, with w,j - w,„j or a left circulant matrix, with w ^ « w,*, (and subscripts 
taken modulo n in both cases). From now on we assume that IK is a left circulant; similai 
results hold in the other case. 

If the y-th slit is open during a particular measurement the total intensity of light passing 
through that slit is 

Vi 

J G{i>)dv (5) 

Since Giv) is a linear function of F(i/) (Eq. (1)) and F{v) is a linear function of the a^ (Eq. 
(4)) this intensity may be written as 

« “I 

*-0 

where the r,-; are constants characterizing the instrument. Thus specifies the amount of 
radiation that is transferred through an open slit at position y, for unit intensity at the k ~lh 
spectral element. We assume that this is only a function of the separation J —k. independent of 
the actual values of j and k. To include more general situations (for example, misaligned 
masks - see Sections 3 and 4) we allow these constants to vary from measurement to measure- 
ment. Thus we assume that if the y-th slit is open during the ;-th measurement, the total 
intensity of light passing through that slit (i.e. the contribution of this slit to the / -th detector 


reading) is given by 


t; “ £ , (7) 

n-O 

where again ihe i „y are constanis depending on the particular instrument. The /-th detector 
reading ti, is given by the sum of all r; for which the corresponding slit is open: 


'I "I 
;~0 

.( -I (I ■'t 

- X I +e, 

\-0 ;-i) 

H-l 

“ X + <-‘i 

k->) 


( 8 ) 


(9) 


where a, is the error in the t >ih reading and the constants c,^ are given by 

II **| 

- X <> -k »‘'-i 
(-1) 

/r-1 

- ^ {since iy IS a leji cirailaiit) 

,-i) 

(i - I 

” 21 {where r ” i-yj 

(-1) 

<i-\ 

c-O 

This is simpler in matrix notation. Let 

V ” (•no.^t 

.1 - (flo.^l . 

e - (c'o.ei, .... 


( 10 ) 


denote column vectors of measurements, unknowns, and errors, respectively, From (9) and 
(10) we obtain 


71 •* Ca + e 
- rH'a + e . 


( 11 ) 

( 12 ) 


where Cls the matrix with (/'.AO-th entry c,^, and 




/f 


. . >»-\ 

r - 

< 1 » 

t.l 

ll • 

• • f'! 


,.i-\ 

[^1 

rr‘ 

I'i • 

•• tr' 


T is called the transfer matrix of the instrument, In an ideal instrument in which distortion, 
aberrations and diffraction were negligible T would be the unit matrix /, In a real instrument, 
however, radiation that should be exiting through a given mask slit often spills over onto neigh- 
boring mask positions, so that diagonal elements of T are reduced in value and otT-diagonal ele- 
ments of the matrix grow. Since radiant energy incident on a detector cannot be negative, 
0 < t,-( < I, Conservation of radiant energy requires that the sum of the matrix elements in 
any row or column of T equals unity, unless dissipative processes such as absorption or 



scattering play a role. If dissipation is significant these sums, 'i-k or 2 

' I, 

than unity and greater than zero, A zero value for any one of these sums would imply a sys* 
tematic blockage in the spectrometer - either by design, or through faulty construction. If no 
light is lost then T satisfies 

TJ « 7*”V - J 

where 7 is a square matrix of ones. 

Equations (ID and (12) are (he basic equations describing the performance of this multiplexing 

spectrometer, They relate the unknowns ao of Fig. 2 to the measurements 

rjQ, , . , via the matrix C - (c^) or via the matrices T and IK Here IF - iw,j) - (iv,^.,) 
is the configuration matrix that describes which slits are open and which are closed during the 
measurements, and the transfer matrix T (Eq. (13)) characterizes the particular instrument but 
is independent of the multiplexing. The entries In T are determined by Eq. (7) and depend on 
the impulse response Hiu) and on the positions 5 q, , . . We shali see below that by suit- 

ably choosing the transfer matrix Tin (12) we can use the same equation to describe the distor- 
tion introduced by a number of different sources of error. 


2.4 Examples 
(1): No Dtllraction 

As a first example we assume a wide entrance slit and no diffraction, so that the impulse 
response is as shown in Fig, 1(a). We take d\ ” 1/2, so that ~ 1 for 

— 1/2 < t'-vo < 1/2, Hiu—vo) “ 0 elsewhere. We also assume that the exit slits have width 
1, so that 0^+1 — 0; “ I and b 1/2 in Eq. (4) and Fig. 2, The slit function 5'(fo;0^,0/+i) is 
given in Fig. 1(c), or analytically by 

. vo-0y4*Y if 0;-j^vo<0,+j 

5(i>o;0y. * “V()-i-0^-fY if ^ *'0 ^ 

0 otherwise 


(14) 

If the y'-th slit is open during, the /-th measurement, the total intensity of light passing through 
that slit is found by substituting and F{i>) into Eq. (2), to give 

rj - r (vq— 0y4*“){(2y_i + (fly — flj_i){j'o"0y+y)lc/i'o 

«y^•3/2 

+ r i.-v^+Oj-\-^)[aj-v{aj^x-aj){.us)—dj—^)]dv(i 
«^*l/2 ^ ^ 


- -rioj-i + 4<7j -F Uj+i) 


(15) 


Therefore the constants tj-n in Eq. (7) are given by 

, 4 , , I 

'0 “ 6 - '-1 ” ■ 
tj =■ 0 otherwise . 

(Note that in this instrument the tj are independent of i.) The transfer matrix Tfor this instru- 
ment is therefore 


4 

1 

0 • 

• 0 

1 

1 

4 

1 • 

> 0 

0 

0 

1 

4 • 

• 0 

0 

6 

6 

6 ■ 

• 4 

i 

.1 

0 

0 • 

• 1 

4 



( 16 ) 
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Exanipte (2): Full Diirractioii 

The second example assumes a very narrow entrance slit, so that the impulse response is 

- 2 swc*(27t{v , (17) 

as in Fig, Kb) with di - !/2, (The constant 2 is chosen so that the total area under 
is one). Again we assume the slits have width I, — 0^ ^ \, i - 1/2, Now the contribu- 
tion of the /-ih sill to the /-th measurement is (from Eqs, (2), (4), (17)): 

’■y' " Z S {«<■ + (a,+i . 

'-0 «,*|-l/2 ^ 

• r 2 s/nc^27r(i^ — 


After some algebra this becomes 

t; - iiiOj + (r'lfly-t + r'-ifly4-i) + • ■ • 

(i.e. Eq. (7)), where the coefilclents are given by 

in 



(18) 


fix) - y(.V + y)^ if -y < ,\f < -y , 

/(.V) - (19) 

fix) - |(.v if I < ,V < I . 

Numerical evaluation of the integrals (18) shows that the transfer matrix T of this instrument is 
a symmetric circulant matrix with first row equal to 

(0.6667, 0.1482, 0.0080, 0.0031, 0.0017, 

0.0010, 0.0007, 0.0005, 0.0004, 0.0003, 

0.0003, 0.0002 0,0017, 0.0031, 

0.0080, 0.1482) 


( 20 ) 


correct to four decimal places. Notice the very slow decay of elements off the main diagonal, 
illustrating the pronounced spreading effect of this impulse response. 

in general whenever H{p) is an even function, i.e, satisfies H{—v) - H{v)y then Tis a 
symmetric circulant — a matrix of the form 


"0 

Ml 

Ml 

... Mj 

M: 

Ml 

»l 

Mo 

M| 

• » » W4 

M 3 

M; 


Ml 

Mq 

• • • Mj 
« « • 

lU 

Mj 

“1 

Ml 

Mj 

■ ■ ■ Wi 

M| 

Mo 


( 21 ) 
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2.6 Optiratlon as a ,Vlotiochromu(or' 

If the same instrument is operated as a conventional monochromator, n measurements are 
made, where in the /*th measurement a single exit dii extends from 0, to The basic 
equation is now 

- Ta + e . (22) 

obtained by setting If' « / in (12). 

2.7 Recovery of Spectrum 

To estimate the spectrum from the measurements we proceed as follows. We assume 
(again this is the ideal case) that e„ the detector error in the /-th measurement, is a random 
variable with mean zero: 

E[e,\ - 0 , (23) 

that the errors in distinct measurements are uncorrelated: 

E{e,ej\ “ 0 , i j , (24) 

and that 

fie,’) >“ cr^ (25) 

where the variance tr^ is independent of the amount of radiation falling on the detector. Under 
-these assumptions if (12) holds then the best estimate of the fl,’s is given by 

a - IK-'r"'T) (26) 

« a -I* IK"' r~'e . (27) 

This is best in the sense of being that linear unbiased estimate which minimizes the average 
mean square error 


,«l''f'£((a,-a,)2) 

« ,to 

- ~ Trace\{TW)-'[{TW)-^)^ 

n 

I 

- — X {sj/m of squares of elenienis of (TH^)”'] (28) 

n 

- see Refs. 3,18. Then e is a measure of the accuracy of the experiment, and should be made 
as small as possible. 

2.8 Computing the Spectrum. Inverse Matrices 

The spectrum is estimated by (26). Usually IK"' is known e.xplicitly,^ but problems may 
arise in finding T~', since in general the inverse of say a 255x255 matrix is both ditficult to 
compute and awkward to store. 

If Tis a circulani matrix, as is often the case, then the inverse is fairly easy to find and to 
store “ see below. If T is Toeplitz (constant along diagonals) but not a circulani. then the 
inverse can still be found but with more difnculiy.''^"^’ But usually if F is not a circulani then 
the best method of computing the spectrum is not to attempt to find F~' but to directly solve 
the system of equations 

T) ■=* riKa (29) 

for a. This is particularly straightforward in the important case when IK"' is known and Fis a 
band matrix, i.e. the only nonzero entries in F are those within a fixed distance of the main 
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diagonal. 

The following are two methods for finding the inverse of an invertible circulant matrix C 

with first row co-^i Method often gives an explicit formula for the inverse: 

Let c(.t) “ Ci)+C|.v+ • • ■ Find the unique polynomial 

<i(x) “ (/o+(/i.v+ • • • such that c(.v)i/(.v) s 1 (mod x" - 1). Then C"‘ is the cir* 

culant with first row da, .. . , J„-|. Method (II) (which may be new) is an efilcient algorithm 
for finding the inverse of a symmetric circulant with a computer: For large n the m-ih entry d,„ 
of the first row of the inverse approaches 

-L J cos m&dO QQj 

^ ® Co + 2 ^ c,cos rQ 
/-I 

In calculating (26) it Is useful to remember that any two right circulants commute, l.e. 
satisfy - BA, any two left circulants commute, and any left circulant commutes with any 
symmetric right circulant. 

Example (1) (cent.) 

The inverse of (16) was found by Method (I). First define the sequence of integers 


Ao “1. Ai - -’4 

b„ “ “4A„.-I ~ b„ -2 , n 2 , (31) 


This is Sequence 1420 in Ref. 23. The first few terms are 

Ao 

I -4 15 -56 209 -730 

An explicit formula Is 

b„ “ |((3-2V3)(-2W3)»+ (3+2V3)(-2->/3)''l 
6 


and also 


A,-_i - A„A„_2 “ I , /I ^ 1 . 

Then the inverse of (16) is given by the symmetric circulant (21) with 


_ 3A„-| 

2A„_| 4- A„-2 + i 

_ 3(A„-,„t -1“ ^.-|) . > j 

2A„-i + A„_2 4- 1 ’ ' 


(32) 

(33) 


(34) 

(35) 


(This may be verified by using Method (I): we omit the rather lengthy details.) For e.xample 
when /I =• 4 

7-2 1-2 

-2 7-2 1 

1-2 7 -2 * 

-2 1-2 7 

For large n the second term in (32) dominates and we can approximate u, by 

l)'V3 (2 -^/3)' 


4 10 1 
14 10 
0 14 1 
10 14 


(36) 
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for 0 ^ «/2, Independent of n. The approximation is very accurate, in fact 

if/f>20. 

Table I gives the values of ii', correct to three decimal places. Note that u\ approaches 0 
rapidly as / increases. Also 


del r = 


2 + 75 
6 


(37) 


Table I 

First Row of Inverse Matrix of (16) 


i 

It', 

0 

1.732 

1 

-0.464 

2 

0.124 

3 

-0.033 

4 

0.009 

5 

-0.002 

6 

0.001 


Therefore for example if the instrument is being operated as a monochromator, so that 
Eq. (22) applies, the best estimate of the fl,'s is given by 

a - . (38) 


and for n ^ 20 we can write Eq. (38) as 

a, “ 2/oT?i + + • • ■ + w'6(tJ,_.4+7J,-4) 

with an error of less than 10"^ where the ii‘j are given by Table 1, and subscripts are read 
modulo n. The average mean square error « for this mode of operation can be found from Eqs. 
(28), (36) and is 

e ~ 2V3 cr* . (39) 


Incidentally, if the correction matrix is' not applied, the result is a coarsening of the spec- 
trum. Consider an input spectrum with a single sharp line, corresponding to 
a - (0,0,0, 1,0,0,...) as shown in Fig, 3(a). Then from Eq. (22), ignoring e, 


7j - ra - (0.0, 1/6, 2/3, 1/6, 0,0,...) . 


If we ignore T~' and estimate a by tj (instead of (38)) the resulting spectrum is shown in Fig. 
3(b). 


Suppose the same instrument is operated as a multiple.xed spectrometer, so that Eq. (12) 
applies. The argument used in the Appendix of Ref. 15 can be modified to show that for any 
(0,l)-matrix IK we have 


^ 10.34 cr- 

e ^ 

n 


n large , 


(40) 


(using Eq. (37)). If IK is an S-matrix'^ then it can be shown that 

13.36 0-- 

, /) large , 


€ = 


/! 


(41) 
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whlch is reiisotiubly close to (40) and may well be the smallest e that can be achieved. Compar- 
ing (39) with (41) we see that 

/n,s,e. mth intiltiplexiiif’ __ 4.00 
m.s.e, wtiliotK nuiltipluxiiis n 

and in fact for this choice of l/lt can be shown that (42) holds for any matrix T. 

E.xamplc (2) (cont) 

The inverse of (20) was found by Method (il), and is a circulant with first row which 
approaches 

(1.6683, -0.3820, 0.0705, -0.0183, 0.0017, 

-0.0017, -0.0006, -0.0006, -0.0004, -0.0003,...) (43) 

correct to four decimal places, for targe /t. Again the elements off the main diagonal decay very 
slowly. Without multiple.xing the average mean square error for this instrument is 

e =!= 2.93 0*^ . 

while multiplexing with an S-matrix gives 

11.74 

€ =s , 

n 

and again (42) holds. 

Of course for any instrument we can ask the following question. The impulse response 
H{v) determines the transfer matrix T, as in Eq. (13). Then which configuration matrix W 
minimizes the average mean square error (28)? In Ref. 15 we studied the performance of vari- 
ous M^s in the very special case when T is the identity matrix (thus Ignoring the spreading 
effect of Hiv)). The more general question for an arbitrary Hiv) or T remains unsolved. 

3. Errors Occurring with a Mask which is Moved in Steps 

This section deals with some classes of errors that can occur with a mask which is stepped 
to the next position between measurements. Most of these can be described by a suitable 
modification of the transfer matrix Fin the basic equation (12). These errors may also arise in 
instruments with a continuously moving mask (see § 4), The resulting Fis then the product of 
the matrices described in §4 and those below. The echo effects produced by slits which are uni- 
formly too narrow (or too wide) were described by Tai et al.*'* 

3.1 Faulty Mask Alignment 

Sometimes the position of an encoding mask is systematically displaced by some small dis- 
tance 5 from its correct position (see Fig. 4). Plankey et al,'* have described the effects pro- 
duced by deliberate misalignment of a spectral encoding mask by 0.5 and 2.5 slit widths. .As 
one would expect, the main result is an apparent spectral shift by 0.5 and 2.5 resolution ele- 
ments, respectively, and an increase in spectral noise. The e.Nact form of this noise may be 
determined as follows. We assume that J5 is small compared to the slit width. Then the contri- 
bution to the ;-th detector reading from light passing through the >-th slit is now given by 

CU>)dv (44) 

instead of Eq. (5). For the instrument of Example (1), this becomes 




- u • 


1/ -U2“A ^ ^ 

f ~5; + -^ + 5 + ;»()){a,„i + (fl/ —fl, -iJCt'o - + ~)l(/i^q 

« -1/7 i ■ i 

«,*\n 

+ r (0y ^ - 5 — V^){a^_| + (a^ — a,_l)(l»g — , + '~))t/;/g 

»».+|/2-4 ^ ' * 

«,+J/2~4 

’ r x 1 

4" J (0y + — “* 3 ~ vq)[oj + (fly+i — Oj — —)\dv,f 

+y 6^(flj_i — 2ay+a^+l) + -I* 5^(fly_j — 3fl; -|+3o^ -a,-i) (45) 

Hence the transfer matrix fls an asymmetric circulant with first row equal to 

4-(4-6S^4-35^(l-5)^0 0,5^ 1+35+35^-33^) , (46) 

g 

Thus a sharp spectral line is dispersed asymmetrically into the adjoining resolution elements. 
Furthermore a sharp line at one end of the spectrum affects elements at the opposite end of the 
spectrum. 

To deal with this problem the mask should first be repositioned so as to minimize the 
misalignment. The remaining error can then be removed by using the correct T"' in Eq. (26), 

3.2 Differences Between Slit Width and Step Size 

When a mask contains large number of slits, a similarly large number of steps must be 
taken in order to move the encoding mask from one extreme position to the other. Over such 
a range of steps, small differences between the step size and the encoding slit width can accu- 
mulate, For example, a typical spectral mask might have a slit width of 0.1 mm and be 255 
slits wide. Over this total width of 2.55 cm, the mask must be stepped with suflicient precision 
that its final position is precisely one step short of cycling into the initial configuration. Thus 
the systematic error in the mask motion has to be less than 4x10“^ mm per step, for the final 
position to be within one-tenth of a slit width from the intended location, 

If such precision cannot be attained then sharp spectral lines will take on a broadened 
appearance. Figure 5 shows a computer simulation of this effect. The analysis of this distortion 
is very similar to that In Sec. 3.1. 

Suppose the initial position of the mask is correct, and the final position is nS behind the 
correct position. Then Eq. (5) must be replaced by 

r;- / G{u)di> , 147) 

The transfer matrix Tfor the instrument of Example (1) has /-ih row equal to 

(0 0,(/A)M+3/A+3(/A)2-3{(A)^ 

4-6(/A)^+3(/■A)^(l-tA)^0 0) , 148) 

with the entry 4 -6(/A)^+3(/A)^ on the main diagonal. 



3.3 Excessive Gap Between Encoding and Blocking Masks 

ir the blocking mask and encoding mask are mounted too dose together, there is a danger 
that one of the opaque dements of the encoding mask may catch the edge of the blocking 
mask, To avoid this there must be a gap between the two musks. If ihis gap is too large, how- 
ever, some of the radiation that strikes the blocking mask obliquely will pass into an encoding 
slit that actually was meant to be obscured by the blocking mask and therefore inaccessible to 
any incident radiation. Then a spectral line imaged on either e.xtreme end of the blocking mask 
will produce a false spectral line at the opposite end of the spectrum. This effect is most pro- 
nounced in fast optical systems (with small focal ratios). 

As long as the encoding (rather than the blocking) mask is placed in the plane where the 
sharpest focus is obtained, this effect is not very serious because it oniy affects the e.xtreme ele- 
ments at the two ends of the spectrum. In contrast the mask misalignment mentioned in §3.2 
produces distortion over the entire spectrum, 

3.4 Nonlinearitics 

Any of the versions of Eq. (12) derived in this paper will only be valid as long as there is 
a linear relation between the input and the detector readings. Eventually this linear relation will 
fail, as the number of spectral elements increases, or as the resolution is increased, for exam- 
ple, Equation (12) must then be replaced by a system of simultaneous nonlinear equations. 
Although numerical techniques for solving systems of nonlinear equations are available,*^' lit- 
tle work has been done so far to determine the e.xaci form of the equations that will be needed 
to replace (12), 

Some nonlinearilies can be compensated for very simply: 

(a) A nonlinear response of the detector to increasing amounts of radiation can be elim- 
inated through judicious calibration which permits nonlinear (compensating) scaling of the 
actual detector output before further data processing is undertaken. 

(b) A nonlinear wavelength or wavenumber response - when Giu), the spectral display 
function along the exit plane, is not a linear function of v - can be compensated for by suitably 
plotting the final spectrum in a way that takes the actual wavelength calibration of the instru- 
ment into account. 

(c) However, if the detector response at one wavenumber is a nonlinear function of the 
other wavenumbers present, compensation becomes very difficult, and the best procedure may 
be to solve for the spectrum iteratively, by successive approximations. 

4. Distortion Introduced by a Continuously Moving Mask 

In this section we analyze the distortion introduced when the slits move continuously 
across the exit focal plane instead of being discretely stepped between measurements. (An 
alternative analysis, which applies also to imagers, has been given by Gottlieb.' However, the 
matrix approach given here seems simpler.) We shall see that in many cases the basic equation 
(12) still holds, with fgiven by (13), where the entries are obtained from Eq. (7). Again 
the spectrum is estimated by a - and the effect of T"' is (theoretically at least) to 

eliminate the distortion caused by the moving slits. 

4.1 Derivation of Basic Equation 

As mentioned in §2.3 we assume that the configuration matrix is a left circulani; 



Wq 

... 

w„-| 


w, 

IVj ... 

Wo 


W2 

Wi • - • 

Wl 


sv,,-! 

Wo ■ • * 

'V-I -2 


In many cases the S-matrices^' can be arranged so as to have this form. An example with 



« •- 7 is given in Eq. (50). 


I I 

I < 


1 

0 

1 

0 

0 


0 

\ 

0 

0 


I I 


0 I 

1 0 

0 0 

0 I 

1 I 

I I 

I 0 


0 0 

0 I 

1 t 

I I 

I 0 

0 I 

1 a 


(50) 


Then instead of using n separate masks of length n -- one for each row of (49) - we can use 
one long mask of length 2n which is moved continuously across the exit plane. Figure 6 shows 
the mask of length 14 corresponding to (50). The extra half-stii at each end of the mask is 
necessary to avoid errors at the ends of the spectrum, if w„_| « 0 in (49) the e,xira half-slits 
are opaque. The mask is periodic with period tu 

Ideally the slits have width lb “ Q,+\ — 9j. Let T.) be the time for each measurement. 
When operated correctly the mask should move at a constant velocity 2bfT(^, so as to move one 
slit width in the time taken to make a measurement. 

In general, as in the previous section we let rj denote the contribution towards the /-th 
detector reading from light passing through an open slit which crosses the y -ih segment 
(see Fig. 7). We may write 


r: 




(51) 


where i/,^, are the endpoints of this slit at the start of the /-th measurement. 

are liie endpoints of the slit at the conclusion of the / -th measurement, and liv) is the 
length of time during which light at wavenumber t/ can pass through the slit, Thus Ki>) can be 
read off Fig, 7; it is the length of the intersection of a vertical tine at position i> with the shaded 
region which is the path swept out by the slit. By combining (S) and (51) we again get the 
basic equation (12), 


4.2 An Ideal Case 

Figure S shows the ideal case, when the mask has the correct dimensions, is properly 
aligned, and moves at the right speed. In this case 
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IM - , 0, ~ b it y 

l(p) » 'J.*t^b-u)ro/'2b , 0,^1 - b K . ■': C,^, + b . (52) 

As an example, suppose the instrument of Example (1) is being operated in this ideal mode. It 
follows from Eqs. (4), (5U, f52) that 

Tj *j|^(a^_2+76a^~i+23Qa_,+76a^4.i+a^4.}) (53) 

Since this is of the same form as Eq. (7) we see that (U) still holds, where the transfer matrix 
TIs now a circuiant with first row equal to 

—^(230,76,1,0,0 0,0,1,76). (54) 

Actually this is an approximation. Since the input spectrum F(v) Is in fact not periodic as we 
have assumed, Eq. (S3) must be modified for i - 0, 1, 2,/i -3,/i ~2,/i -I, and the top right- 
hand and bottom left-hand corners of Tare missing. Thus 7'is a Toepiitz*'^ matrix, not a circu- 
iant, However, it seems worthwhile changing f to a circuiant and accepting the resulting distor- 
tion in the ends of the spectrum, in order to obtain a matrix which has a manageable inverse. 
The inverse of the circuiant (54) Is easily obtained, e.g. by Method (II) of §2,3, The inverse is 
a circuiant with first row which approaches 

(2.213, -0.826. 0.299. -0.108, 0.039, 

-0.014, 0.006, -0.002, 0.001, 0, 0, 0....) (55) 

correct to three decimal places, as n increases. For /» "• 30 (55) is already valid to this order of 
accuracy. 

Again the eiTcci of ignoring T is to broaden the spectrum, but now the broadening is 
more pronounced than in Fig. 3. 

4.3 If the Mask Velocity or Silt Width is Wrong 

It is not difficult to modify the above analysis to determine the distortion introduced if the 
mask Is misaligned or if the mask velocity or slit widths are incorrect. It can be shown that, as 
long as each slit has the same (possibly Incorrect) width, Eqs. (12) and (13) still hold. E.g, for 
the instrument of E,xample (1) operating with a misaligned mask moving with incorrect velocity 
we find that 

t/ " fiUj-2 + 4- fQflj 4- (56) 

for certain constants t, (again with suitable modifications for the first and last few measure- 
ments). 

It does not seem worthwhile giving further examples, since each instrument will have its 
own transfer matrix fto be used in Eq, (12), and T should be determined when the instrument 
Is calibrated. 

4.4 Rotating Two-Dimensional Masks 

The above efiects also occur in imagers,^' when the encoding mask may be a two- 
dime-nsiona! mask tnounted on a rotating disk. For example, a slight radial eccentricity will 
cause an efiective radial motion of the mask. The result of this run-out should be a superposi- 
tion of three efiects. There should be a widening of the image along the radial direction, simi- 
lar to the effect discussed in §3.2 above. The overall appearance of the image should become 
noisier, because of the continuous motion, as in §§4,1, 4.2. Thirdly, there should be a single 
cycle of roughly sinusoidal change in intensity along the unfolded two-dimensional chain of ele- 
ments, somewhat similar to an effect described in §5 below. This is because the eccentric 
motion of the mask will move an opaque portion of the wheel into the open frame in the block- 
ing mask, reducing the intensity during that portion of the cycle. 
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4.5 Imugint; a Movinii Source 

If one of ihe sources being observed moves slightly during the measurements, the results 
will be similar to those described in §3.2 and §4.4, depending on whether the motion is along 
the direction of the mask's motion or across it. By and large, the result of a small displace* 
ment, of the order say of one spatial element during the frame time (or complete cycling time) 
of the mask will be to stretch the image along the direction of motion, and to add noise to the 
entire scene. 

5. Effect of Drift in Background Level 

if the intensity of background radiation incident on a spectrometer or imaging system 
varies during a spectral run, the derived spectrum will necessarily be affected, Suppose for 
example that the infrared spectrum of a star is to be observed. If the foreground atmospheric 
emission drifts during the course of the spectral run the instrument will record a corresponding 
drift, and the final spectrum obtained for the star will contain a component which can be 
directly attributed to the drift. Because the spectrum is estimated by a linear operation (Eq. 
(26)), the drift component is simply added to the true stellar spectrum. In this section we 
analyte the magnitude of this effect for several different types of drift, assuming for simplicity 
that T - L (The analysis in the general case is much more complicated and we do not go into 
it here.) 

Let 7) denote the vector of measurements with no drift present (Eq. (12)), let dj be the 
drift in the y*th measurement, 

d - {da,d \, . . . ,fl',_|)^, 

and let 

Tj' •” T) 4* d (57) 

be the actual measurements. The estimate of the spectrum is 

a - fK-'vj' 

+ iy-‘d , (58) 

and we wish to analyze the drift component 

D - IK-'d . (59) 

For concreteness we take to be a symmetric, left circulant S-matrix,^**^ with 

tf/-f „ 5-1 - (25-y) . (60) 

/)+l 

Case CO ” Constant Offset 
Suppose 

d - 61 (61) 

where I - (1, 1,..., 1) Then 

D - 1 . (02) 

Case (il); Sinuosoidal Daft 
Suppose 

df ^ zos 1-n a ij —<b) , y ” 0, 1 /i— 1, (63) 

where a,<b are constants, <t> being a phase shift. The ^‘-th component of the drift is 

fl-i 

' ” 21 ^jOOS 2-ffaij —k —tb) 

7-0 


( 64 ) 



where 


is the first row of Simplifying (64) we obtain 

"* <•/ cos 2:ra(^+ip — /?) , 

where the amplitude Is given by 

= i 

n 


( 63 ) 

(6t>) 


(67) 


for large /i, since (65) is a pseudo-random sequence, and the phase shift /3 Is given by 

-I -I 

22 f;Sin 2TTaj 

tan(27Tfl^) - *;~J . (68) 

22 f;Cos lira, 

;-0 


We conclude that a sinusoidal drift voltage adds to the true spectrum a sinusoid of the same 

frequency as the drift, with amplitude multiplied by 2/'//T, and with a phase shift given by (68). 

Notice that a sinusoidal drift affects the spectrum more strongly than a constant offset of equal 

amplitude. 

» 

Case (ill); A Single Noise Spike 

Suppose dff - 5^ c/y - 0 for y A: Then [Dj] - is a pseudo-random sequence of 

r.m.s. value 


-so; 

" 7 


in 


= 28//I 


(69) 


Thus a single noise spike of amplitude S produces fairly random noise over the whole spectrum 
with r.m.s, value 25/ n. 

Case (iv) 

Conversely, if dj « D has a single spike of amplitude 5 in the A'-ih component. 

These results have been collected in Table II, where the different drifts have been nor- 
malized so as to have unit amplitude. The last line of the table shows the improvement in 
signal-to-noise ratio obtained by considering random noise, obtained from Eq. (42), It is 
interesting that this is approximately equal to the improvement obtained when the noise is 
coherent. Only a constant offset In the data gives rise to a considerably smaller mean square 
error in the final spectrum. 

Similar results are obtained in Fourier spectroscopy. A noise spike in the data will again 
produce a wide distribution of error signals in the final spectrum, A sinusoidal noise com- 
ponent, on the other hand, will produce a single spectral line at a frequency corresponding to 
that of the noise, (The analogous noise in the Hadamard instrument is that considered in Case 
(iv).) 


6. Singular Designs 
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Table II 

Noise in Spectrum Produced by Drift of 
Unit Amplitude in the Data 


Type 

of 

drift 

R.m.s. 
drift 
in data 

Noise 

in 

spectrum 

R.m.s, 

spectrum 

noise 

Ratio of 
r.m.s. spectrum 
noise to 
r.m.s. drift 

Constant 

oll'sct 

1 

Constant 

ofiset 

2/ a 

Vn 

Sinusoid 

mwm 

Sinusoid 

v'2//j 

IN n 

Single 

Spike 

I/^/^ 

Eq. (69) 

2//1 


Random 

• 

Random 

- 

2Nn 


6.1 More Measurements Than Unknowns 

It is sometimes desirable to design the e.xperiment so that the nui:<ber of measurements, 
p, exceeds the number of unknowns, //. The purpose in doing this is to ensure that in case 
some measurements are lost (for example if a cloud passes over an observatory during an astro- 
nomical observation) the spectrum can still be calculated.’** 

The basic equation (U) relating the vector of measurements tj - ( 110.^11 i}p_i)*’to 

the vector of unknowns a becomes 

7} - riPa + e . (70) 

where IP" Is a pXn (0,1) configuration matrix describing the experiment, and fis a pxp transfer 
matrix characterizing the instrument. The iij) -ih entry of Fis 0 < /,y /? -1 (see Eq. 

(7)). 

The best estimate for a is now-*-*** 


a - (FIK)> , 

(71) 

where (FIK)'^ is the generalized inverse, given in this case by 


(FIK)-^ - (fK*'F*'FIK)">(FIK)*‘. 

(72) 

The average mean square error is then^ 


« - ~ Fwel(IK^F*'FH')'') . 

(73) 


« 


A mask of this type suitable for use in a spectrometer can be obtained by taking If' to be the 
first /I columns of a /jxp circulant S-mutrix. This can be accomplished by means of a blocking 
mask which only exposes a frame of tt slits at a time. 

6,2 More Unknowns Than Measurements 

Suppose an experiment has been designed to make n measurements in order to determine 
a spectrum with /( unknown components, via Eqs. (12), (26), but is prematurely terminated 
after only p < n measurements have been made. In some cases it is still possible to say some- 
thing about the unknown spectrum. For example, suppose IK - S,, and F - /. Then tj «= 5,, a, 
and the sum of all n measurements is equal to (/j4-l)/2 limes the sum of the unknown spectral 
components! 

i - i(„+l) i a. . (741 

i-\ (-I 

Suppose the last measurement i,, is missing: what can be said about its possible value? Since 
is the sum of those a, for which there is a I in the last row of we have 














( 75 ) 


T?„ < 2 • 

.-1 

From (74) and (75) wc obiain ihe bounds 

0 < ^ Z ’l; • ^76) 

A similar result can bo given if two measurements are missing. In some cases the lower bound 
in (76) can be improved by considering the dlircrences ~ V, (dr suitably chosen / and J, 

An analysis of this type can also be given for other multiplexing schemes, including 
Fourier spectroscopy. A single missing measurement never entails complete loss of informa- 
tion about the spectrum. The uncertainty however increases rapidly as the number of missing 
data points increases. 

6.3 Correction Procedures 

A similar situation arises when one or more measurements are lost because of a burst of 
noise. A single noise spike of amplitude 5 produces random noise over the whole spectrum 
with r.m.s. value 25/«, as we saw in §5. 

A large noise spike can often be recognized by examining the other measurements 
obtained in the same run. When the spectrum is continuous, or contains a large number of 
intense lines, the individual data values do not greatly deviate from each other. Only if one or 
two spectral lines dominate do wide deviations occur, For a polychromatic spectrum, then, 
large noise spikes tend to .appear as well defined, extraordinary points. 

In the laboratory we have tended to remove these spikes and replace them by the mean 
value of adjacent data points on either side. The theoretical justification for this procedure is 
somewhat questionable, but wc find in practice that the spectra obtained by means of such 
corrections are rather good approximations to the' expected forms. The procedure is illustrated 
in Figs. 9-12, Figure 9 shows a spectrum obtained in the usual way with a 255x255 S-matrix. 
Figure 10 shows the elTcct of a noise spike in the last mc.isurement, while Fig. 11 shows the 
result of replacing this and four adjacent data points by a straight line joining the 250-th and l- 
st data points. Similarly in Fig. 12 the last 15 measurements have been replaced by 15 points 
linearly Interpolated between the 240-ih and 1-st measurements. 

A comparison of Figs. 9 and 1 1 suggests that linear interpolation between values spanning 
the gap of missing data entries produces relatively little distortion. A somewhat better pro- 
cedure might be to use a quadratic or higher degree curve for the interpolation. 

Figure 13 uses the same data as in Fig. 9, except that a correction has been inserted to 
remove the negative echo caused by faulty slit deposition (or etching), as in Tai et al. (1975). 
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F!i;uru Captions 

Fig. I, Impulse response Hiu—vn) (a) when dilTraction can be ignored, and (b) when 
dilTraction dominates, (c) and (d) are the corresponding slit functions. tSee Examples tl) and 
(2) in §2.4.) 

Fig. 2. A piecewise linear input spectrum f (i^). a, is the intensity at the midpoint of the inter- 
val 

Fig, 3. Showing how a sharp input (a) is coarsened (,b) if the correction matri.x T*”' is not 
used. 

Fig. 4. Faulty mask alignment. 

Fig. 5. Computer simulation of the errors produced by a systematic dilTcrcnce between slit 
width and step size. The final position of the mask after 255 steps is displaced by one slit width 
from the correct position. An input spectrum with a single sharp spectral line has been 
broadened to a width of two spectral elements, and noise has appeared across the whole spec- 
trum. 

Fig. 6. Mask designed to be moved continuously. This is a mask of length 14 corresponding to 
the lFmatri.x of Eq. (50). Elements 1-6 are repeated as elements 8-13 respectively. The e.xtra 
half slit at each end, corresponding to clement 7, is needed to avoid errors at the ends of the 
spectrum. 

Fig. 7. Path swept out by an open slit crossing the J -th segment () during the t -th detector 
reading. Time runs downward in this diagram, and the mask moves to lhe^ right. 

Fig. 8. The ideal case when the mask has the correct dimensions, is properly aligned, and 
moves at the right speed. .\gain lime runs downward and the mask moves to the right. 

Fig. 9. Uncorrccted 255-element spectrum of a laboratory source showing the 1.7/j, mercury 
vapor lines. There has been no data processing other than applying the inverse Hadamard 
transformation. 

Fig. 10, Same as Fig. 9 e.xcept that the last data point, of height l.lb.xlO^ has been replaced 
by a noise spike of height 4x10^, The main spectral features remain, but noise is added to the 
entire spectrum. 

Fig. 11. Same as Fig. 9, e.xcepl that we have simulated losing the last five data points and 
replacing them by a straight line joining the 250-th and 1-st data points. 

Fig. 12. Same as Fig. 11, but now the last 15 measurements have been replaced by 15 points 
linearly interpolated between the 240-th and 1-si data points. 

Fig. 13, The data used here is the same as in Fig. 9, e.xcept that a correction has been inserted 
to remove the negative echo caused by faulty slit deposition (or etching), as in Tai et al. 
(1975). 
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